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ABSTRACT 

We present the first investigation of the Blue Straggler star (BSS) population in two isolated dwarf 
spheroidal galaxies of the Local Group, Cetus and Tucana. Deep HST/ACS photometry allowed 
us to identify samples of 940 and 1214 candidates, respectively. The analysis of the star formation 
histories of the two galaxies suggests that both host a population of BSSs. Specifically, if the BSS 
candidates are interpreted as young main sequence stars, they do not conform to their galaxy's age- 
metallicity relationship. The analysis of the luminosity function and the radial distributions support 
this conclusion, and suggest a non-collisional mechanism for the BSS formation, from the evolution 
of primordial binaries. This scenario is also supported by the results of new dynamical simulations 
presented here. Both galaxies coincide with the relationship between the BSS frequency and the 
absolute visual magnitude My found by iMomanv et a l. (2007). If this relationship is confirmed by 
larger sample, then it could be a valuable tool to discriminate between the presence of BSSs and 
galaxies hosting truly young populations. 

Subject headings: (stars:) blue stragglers Local Group galaxies: individual (Cetus dSph, Tucana 
dSph) galaxies: evolution Galaxy: stellar content 



1. INTRODUCTION 

Blue Straggler stars (BSSs) were first identified in th e 
Galactic globular cluster (GGC) M3 bv lSandaed (fl953h . 
as a group of stars bluer and brighter than the cluster 
turn-off stars, thus being more massive than the stars 
currently evolving at the turn-off (TO). Since then, BSSs 
have been identified in a variety of stellar systems in 
very different environm ental conditions including: th e 
field (jCarnev et a l.||2005| ), open clusters fLiu et al.l l2008h. 
globular clusters (e.g.. | Piotto et ail 120041), an d dwarf 
galaxies (|Lee et al.l 120031: IMomanv et al.l 12007ft . Their 
puzzling position in the color-magnitude diagram (CMD) 

1 Based on observations made with the NASA/ESA Hubble 
Space Telescope, obtained at the Space Telescope Science In- 
stitute, which is operated by the Association of Universities for 
Research in Astronomy, Inc., under NASA contract NAS5-26555. 
These observations are associated with program 10505. 

2 Instituto de Astrofi'sica de Canarias, La Laguna, Tener- 
ife, Spain; monelli@iac.es, carme@iac.es, antapaj@iac.es, shi- 
dalgo@iac.es. 

3 Dcpartamcnto de Astroffsica, Universidad de La Laguna, 
Tenerife, Spain 

4 INAF-Osservatorio Astronomico di Collurania, Teramo, 
Italy; cassisi@oa-teramo.inaf.it. 

5 Univcrsita di Milano Bicocca, Dipartimento di Fisica 
G.Occhialini, Piazza della Scienza 3, 1-20126, Milano, Italy; 
mapelli@mib.infn.it 

"Institute for Astronomy, University of Edinburgh, Royal 
Observatory, Blackford Hill, Edinburgh EH9 3HJ, UK; 
ejb@roe.ac.uk 

7 Department of Astronomy, University of Minnesota, Min- 
neapolis, USA; skillman@astro.umn.edu. 

8 Dominion Astrophysical Observatory, Herzberg Institute of 
Astrophysics, National Research Council, Victoria, Canada; 
peter. stetson@nrc-cnrc.gc.ca. 

9 Department of Physics, Institut fur Astronomie, ETH 
Zurich, Zurich, Switzerland; lucio@phys.ethz.ch. 

10 Institut fur Theoretische Physik, University of Zurich, 
Zurich, Switzerland; lucio@physik.unizh.ch. 

11 Kapteyn Astronomical Institute, University of Groningen, 
Groningen, Netherlands; etolstoy@astro.rug.nl. 



suggested that they do not fit in the traditional scheme 
of stellar evolution for single stars. In fact, in order to 
populate that region of the CMD, these exotic objects 
must have experienced a physical process able to allow 
them to still stay in the core H-burning stage despite of 
their mass and the clust er age (for a detailed historical 
review, see IStrvkerlll993l) . 

Presently, two different physical mechanisms are fa- 
vored for explaining the formation of BSSs: i) co- 
alescence throug h direct stellar collisions (COL-BS, 
iHills fc Davlll976ft . and ii) evolution of primordial binary 
systems, in which mass transfer between the two compo- 
nents allows th e rejuvenation of the secondary compo- 
nent (MT-BSS. lMcCrealll964[ ). It is also widely accepted 
that these two scenarios are not in competition, but they 
could both occur simultaneously in a stellar system. In 
fact, the coexistence of BSSs formed through both the 
collisional and mass transfer channels is commonly in- 
voked to interpret the radial distribution of BSSs in many 
GGCs. The higher central concentration is explained by 
an efficient rate of stellar collisions in the densest regions, 
while the rise at larg e distance would be due to MT-BSSs 
(jMapelli et al.ll2006t iFerraro et al.ll2009f ) surviving in the 
outskirts. 

While the scenario for BSSs in Galactic clusters is well 
settled, at least from the observational point of view, in 
contrast, this is not true for dwarf galaxies. In fact, even 
though the first indication for the presence of BSS can- 
didate s in the Sextans dwarf dates back to iMateo et al.1 
(1991), and, since then, several observational findings 
suggest t he presence of BSS candidates in several dSphs 
(see, e.g- JMapelh et al.1l2009l and references therein), we 
are still faced with the problem of understanding if these 
candidates are genuine BSSs or normal core H-burning 
stars belonging to an, albeit sparse, intermediate-age 
stellar population. 
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TABLE 1 

The distance modulus, reddening, 
and candidate bss properties. 



Cetus 



Tucana 



(m - M) 


24.46 ±0.12 


24.74 ±0.12 


E(B-V) 
N flS a 


0.028 


0.031 


940 


1214 




90 


123 


Fbs c 


0.05 ± 0.05 


0.02 ± 0.04 


Kc 


3 


6 



a Number of candidate BSSs 
k Number of evolved BSSs in the core He- 
burning stage 

c Logarithmic frequency of BSSs with re- 
spect HB of stars 

Number of Anomalous Cephcids 



Recently, a few studies have focused on the proper- 
ties of BSSs in nearb y dSp h galaxies in the Local Group 
(LG). lMapelli et all (|2007l) presented wide field data for 
the galaxies Draco and Ursa Minor, concluding that the 
population of BSSs is, i n both systems , comp atible with 
the MT-BSS scenario. iMapelli et all (|2009T ) discussed 
the properties of BSSs in Sculptor and Fornax. Sculptor 
presents similarities with both Draco and Ursa Minor, 
as its BSS population is also compatible with a MT-BSS 
population. On the contrary, the relatively strong cen- 
tral radial concentration of the BSS candidate stars in 
Fornax favors the presence of an intermediate-age stel- 
lar population (in agreement with previous studies of the 
star formation history (SFH) of Fornax). 

Unraveling the nature of the BSS candidates has im- 
portant implications for our understanding of the SFHs 
of these galaxies. For example, whenever it is possible to 
discriminate between a sequence of BSSs and a young MS 
population, it is possible to set a constraint on the age 
of the last star formation event in a galaxy. Therefore, 
the BSS populations in Draco, Ursa Minor, and Sculptor 
suggest that no star formation occurred in the last few 
Gyr in these g alaxies. 

Intriguing ly, iMomanv etHI (|2007| ) found a statisti- 
cally significant anti-correlation between the frequency 
of BSSs and the absolute visual magnitude (My) of 8 
nearby dSphs. The anti-correlation is valid over a large 
range of luminosities, and p resents a different slop e than 
that of GGCs. In addition, IMomanv et al.l (|2007f ) ruled 
out the possibility that collisional binaries could have 
contributed to the observed BSS populations and con- 
cluded that the BSSs in dwarf spheroidal galaxies are 
mainly formed via mass exchange in primordial binaries. 

In view of the important implications of the presence 
of a sizable sample of BSSs in a dwarf gal axy, it is worth- 
while to apply the approach outlined by Moman v et al.1 
(2007) to other dwarfs. The hope is to eventually pro- 
duce independent constraints useful for addressing the 
true nature of the BSS candidates. 

In this paper, we investigate these issues further by an- 
alyzing the properties of candidate BSSs in the isolated 
dSphs Cetus and Tucana. The data sets are part of the 
LCID project (Local Cosmology from Isolated Dwarfs), 
aimed at recovering the full SFHs of six isolated LG 
galaxies: Cetus and Tucana (dSph), LGS 3 and Phoenix 
(dlrr/dSph), IC 1613 and Leo A (dlrr). The plan of the 
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Fig. 1. — The observed CMDs for the two dwarfs: Cetus (left 
panel) and Tucana (right panel). The delimited areas represent 
the regions of the CMDs selected for the present analysis (see the 
text for more details). 

paper is as follows: in the next section we briefly describe 
the photometric data set and outline the criteria for se- 
lecting the candidate BSSs in the observed CMDs. In 
§3 we investigate the properties of the BSS populations. 
§4 presents the results of dynamical simulations. The 
discussion and final remarks close the paper. 

2. THE OBSERVATIONS AND THE SELECTION 
OF SAMPLES 

The observations used in this paper have already 
been presented in | Monelli et al.1 (|2010bl Cetus) and 
iMonelli et all (|2010d, Tucana). Here we summarize the 
main points relevant to the following discussion. The 
images were collected with the ACS camera aboard the 
HST, under the project The onset of star formation 
in the universe: constraints from nearby isolated dwarf 
galaxies (PID 10505, PI C. Gallart). A total of 25 and 32 
orbits have been devoted to Cetus and Tucana, respec- 
tively, using the FA75W and F814IU filter pass bands. 
The photometric reduction has been performed usi ng the 
DAOPHOT / ALLFRAME package (jStetsonl 11994). The 
final catalog has been calibrated to the V EGAMAG sys- 
tem u sing the transformations presented in lSirianni et al.1 
(2005). The adopted distance modulus and reddening 
are summarized in Table [TJ together with other quanti- 
ties calculated in this work. 

Figure [1] shows the CMDs of the two galaxies, with 
the regions adopted to select the various samples stud- 
ied superimposed. Note that we adopted the same 
boxes for both galaxies, applying a shift to take into 
account the differences in distance modulus and redden- 
ing. The CMDs of both galaxies are w 1.5 mag deeper 
than the TO, reaching m,F8i4W ~ 28.8 mag. The blue 
plume of candidate BSSs emerges clearly in the CMD, 
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for m,F475W — rnF8UW < 0.5 mag, 26 < tufsuw < 28. 
We defined a box in this region to include the bulk of 
these stars. Note that there are a number of objects be- 
tween the BSSs region and the red giant branch (RGB). 
We decided not to include these stars for two reasons: i) 
it is possible that blended objects from the most popu- 
lated TO region are polluting this part of the CMD; ii) in 
this color range, we expect to find stars evolved off the 
main sequence. Also, we did not include the brightest 
objects (mpsuw ~ 26.2) because it is possible that they 
are extreme horizontal branch (HB) stars. Also note that 
we estimated a negligible contamination from both fore- 
ground Galactic stars (due both to the small area cov- 
ered and the high Galactic latitude of the two galaxies) 
and background galaxi es (thanks to the ca reful clean- 
ing of the catalogs, see iMonelli et aDl2010bf ). Figured] 
also shows the regions adopted to select stars in the HB, 
the RGB, and the progeny of BSSs during the central 
helium burning phase. This last box has been placed 
above the red HB, paying attention not to include the 
AGB clump, which clearly appears in the CMDs of both 
galaxies (jufsiaw ~ 23.4 mag). 

3. ANALYSIS OF THE CANDIDATE BSS 

In the box covering the blue plume, we have iden- 
tified 940 and 1214 objects in Cetus and Tucana, re- 
spectively, which we consider as BSS candidates. Note 
that the present data cover a significant fraction of the 
body of Tucana (tidal radius r t = 3.45'), which was 
centered in the ACS field (field-of-view 3.4'). How- 
ever, Cetus occupies a larger area in the sky (r t ~ 32', 
iMcConnachie fc Irwinl [200Q, and tne sampled field 
covers a smaller fraction. Moreover, the observed field 
was offset from the center so t hat the innermost re gions 
of Cetus were not observed (|Monelli et al.l l2010b[ ). If 
we take this i nto a ccount, following the reasoning of 
iBernard et al.l (|2009f) to estimate the total number of RR 
Lyrae stars, the total number of BSS candidates in Cetus 
is estimated to be ~ 5,000 stars. 

The most important question is to assess the nature of 
these stars. BSSs have been found in very different envi- 
ronments, from the field, to low-density stellar systems 
such as open clusters and dwarf spheroidals, and in most 
globular clusters. They arc therefore common products, 
and it is reasonable to expect to find them in Tucana 
and Cetus as well. Therefore, the basic question we are 
trying to address is: are these candidates genuine BSSs, 
or is there also a significant component of truly young 
stars? In the following sections we will try to address 
this using different observables. 

3.1. Clues from the detailed SFH analysis 

The main objective of the LCID project was to de- 
rive accurate SFHs for a sample of six isolated galaxies 
in the Local Group, in order to understand their evolu- 
tion in a cosmological context. To do this, we adopted 
well-established techniques based on the comparison of 
the observed CMD with a synthetic one. We devoted 
significant effort to investigating the error budget and 
also possible systematics affecting our solution, due to 

12 New preliminary estimates based on wide-field Subaru data 
confirm values rt > 15' (Bernard et al, 2012, in prep.), therefore 
significantly larger than the ACS field-of-view. 



the use of different photometric packages, stellar evolu- 
tions libraries, SFH c odes. The main co nclusion (see 
IMonelli et aHl2010bllcl ; iHidalgo et all 120111) of this anal- 
ysis was that we can rely, within the error bars, on both 
the uniqueness and stability of the derived solutions. 
Therefore, the most important features recovered, such 
as the epoch and the duration of the main star forma- 
tion events, including those at the oldest epoch, and the 
age-metallicity relations are solid results. 

The SFHs presen ted in IMonelli et all (|2010bf ) and 
IMonelli et al.l ()2010c| ) , for Cetus and Tucana respectively, 
show that both galaxies are made of old and metal- 
poor populations. In particular, both galaxies formed 
90% of their stars at epochs older than 9 Gyr ago, and 
the bulk of the stars have metallicities in the range 
0.0001 < Z < 0.001. In addition, both galaxies present 
a well defined age-metallicity relation in the sense of in- 
creasing metallicity with decreasing age. 

The SFH reconstructions indicate a small population 
(< 3% in mass) of stars significantly younger than the 
majority of stars (2 < t < 5 Gyr) but with metallici- 
ties consistent with the oldest stars (Z < 0.0006). Thus, 
these stars do not follow the general age-metallicity re- 
lations in these two galaxies. The analysis of the best-fit 
CMDs show that they nicely populate the blue plume 
region of the CMD, brighter and bluer than the oldest 
MS turn-off. The possibility that the metallicity of these 
stars has been understimated by the code is excluded 
with high confidence level. In fact, higher metallicity 
and an age of « 4 Gyr, would produce also a redder 
RGB sequence that is not observed. On the basis of the 
retrieved SFH, IMonelli et al.l (|2010bH made the hypothe- 
sis that the sequence of blue stars, brighter than the old 
MS Turn-Off, are consistent with a population of BSSs, 
characterized by the same metal content as the metal- 
poor dominant stellar component. 

An alternative interpretation would be that these are 
truly young metal-poor MS stars, formed as a conse- 
quence of gas accretion. However, it seems very unlikely 
that two galaxies, that spent most of their life in isola- 
tion, accreted gas and formed stars of the same metal- 
licity and at the same epoch, in one event that created 
roughly the same amount of mass (3% of the total stellar 
mass). 

3.2. Comparison with stellar evolutionary models 

Figure [2] shows the comparison of the Tucana CMD 
with theoretical tracks f rom the BaSTI databa sJ^I of 
stellar evolution models ([Pietrinferni et al]|2004D . Here 
we note that the analysis of the SFH for the LCID galax- 
ies has been performed by using the same evolutionary 
framework; so the present analysis is fully consistent with 
those concerning the SFHs. The selected tracks span the 
range of metallicities suggested by the SFHs, Z = 0.0001, 
0.0003, 0.0006, and 0.001 as an upper limit. 

The figure discloses some interesting clues. First, the 
black lines in the left panel show the Zero Age Main Se- 
quence (ZAMS) for two metallicities: Z=0.0001 (solid 
line) and Z = 0.001 (dashed). The comparison with the 
observed CMDs shows that both ZAMSs border the blue 
envelope of the BSS sequence, independent of the as- 
sumed metal content. Moreover, the same panel shows 

13 \http://www. oa-teramo.in af. it/BASTI 
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Fig. 2. — The CMD of Tucana with superimposed a selected sample of evolutionary stellar models for the labeled assumptions about the 
initial stellar mass and chemical compositions. Note that, to avoid confusion, the tracks are plotted until the tip of the RGB, and that the 
filled symbols mark the Helium ignition. We have adopted the values of distance modulus and reddening, from Bernard ct al. (2008), as 
summarized in Table [T] Left - The 1.5 and 0.8Mq bracket the entire sequence of candidate BSSs. The two black lines show the ZAMS 
for metallicity z=0.0001 (solid) and Z=0.001 (dashed), which outline the blue edge of the sequence. Central - We show the ranking with 
metallicity of stars of 1.0 (green) and 1.2 (red) Mq. The shaded area encloses the central hydrogen burning phase, until the exhaustion 
marked by the asterisks. Right - Comparison with models with fixed metallicity (Z=0.001) and three different masses (1.0, 1.2, 1.5 Mq). 



that the maximum mass fitting the brightest objects is 
equal to about 1.5 Mq (red lines), again with a negli- 
gible dependence on the assumed metallicity (Z=0.0001, 
0.001 for the red solid and dashed lines, respectively). 
The same figure also shows the location of two theoreti- 
cal tracks corresponding to 0.8M Q stellar models (green 
lines), taken as representative of the stellar structures 
currently evolving at the TO and sub-giant branch of 
both galaxies. The filled symbols show the first point 
corresponding to the central helium burning phase. It is 
evident - as also fully supported by the SFH analysis - 
that for the TO and the SGB, a range of metallicities in 
the adopted evolutionary framework is required to prop- 
erly interpret the observed morphology of the CMD. This 
comparison reveals that the maximum mass of the stars 
in the blue plume is at most ~ 2 times larger than the 
mass of the stars at the TO. This occurrence is in good 
agreement with the properties of BSSs in other stellar 
systems such as GGCs (see, i.e., iGilliland et al.l Il998l 



and references therein). 

The central panel shows the comparison with stel- 
lar tracks for different masses (M = 1.0, 1.2M Q , green 
and red lines, respectively), and different metallicities 
(Z = 0.0001-solid, Z = 0.0003-dashed, Z = 0.0006- 
dot-dashed). The shaded area encloses the region corre- 
sponding to the central hydrogen burning phase, whose 
end is marked by the asterisks. The bulk of BSS can- 
didates clearly occupies this region of the CMD The 
plot clearly shows that a spread in mass is needed to fit 
the observed distribution of blue plume stars. A spread 
in metallicity is also possible but, as shown in the right 
panel of the same figure, values larger than Z > 0.001 
seem unlikely. The 1.0 Mq track is too faint and too red 

14 Although the adopted stellar models are obtained in the 
canonical evolutionary framework, i.e., they do not account for 
the physical processes at the basis of BSS formation, we think that 
in order to have a rough estimate of the metallicity range spanned 
by the BSS candidates, they are still suitable. 
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Fig. 3. — Similar to Figure. [2] but we present a zoom of the CMD 
of Cetus (left panel) and Tucana (right panel), corresponding to 
the core He-burning stage. In each panel the box represents the 
region used to select the candidate progeny of the blue plume stars. 
Four tracks are overplotted, for different mass and metallicity as- 
sumptions. The filled stars mark the locat ion of the anomalou s 
Cepheid variables detected in these galaxies (Bernard ct al. 2009). 



to represent the bulk of the blue plume objects. The 
1.2 Mq still gives a good representation, and can be 
considered as an upper limit both in term in of mass 
and metallicity. In fact, despite the 1.5 Mq still lim- 
its the brightest blue plume stars, the age of this model 
at the turn-off (~1.3Gyr) is in contrast with the find- 
ing of the SFH. On the other hand, higher metallicities 
would shift the tracks to redder color and fainter magni- 
tude, significantly worsening the agreement in particular 
in the RGB phase. Thus, relatively metal-poor and rela- 
tively massive (> 1 Mq) stellar structures are preferred 
to better represent the peak of the distribution. There- 
fore, the metallicity derived for the BSS candidates is in 
good agreement with that estimated for the bulk of the 
older populations, as discussed in i j3.1[ and not in agree- 
ment with the higher metallicities found at the end of the 
initial episodes of star formation in Cetus and Tucana. 
This result can be considered a circumstantial evidence 
supporting the idea that - at least a large fraction of - 
the stars located along the blue plume could be genuine 
BSSs; in fact if they were true young stars formed in a 
late star formation burst, they would have a metallic- 
ity lower than that one could expect of the basis of a 
standard age-metallicity relation. 

Additional information can be provided by the study 
of the BSS progeny in the core He-burning stage. Evo- 
lutionary models for MS stars predict that stars more 
massive than about IMq are located in the red side of 
the H B at a magnitude level depending on their He core 
mass (Castcll ani fc DeglTnnocentil ll995). In particular, 
the right-hand panel of Figure[2]suggests that the evolved 



stars in this mass range all clump in the same region of 
the CMD. This implies that the progeny of BSSs of sim- 
ilar mass, independen t of the formation mechanism (see, 
e.g,. ISills et aL|[2009t) . should occupy a similar region in 
the CMD, redder and brighter than their low- mass coun- 
terparts. This means that, on the one hand, the distri- 
bution of BSS progeny in the CMD provides very little 
information about the mass distribution of the BSS pop- 
ulation above IMq. On the other hand, the ratio of the 
number of stars can be compared with theoretical expec- 
tations based on the evolutionary lifetimes in the core H- 
and He-burning stages. 

Figure [3] shows an expanded view of the HB region of 
the Cetus (left panel) and Tucana (right panel) CMDs. 
Superimposed to both plots are the evolutionary tracks, 
starting with the onset of the central helium burning, 
for two masses (1.0 -blue line- and 1.3 -red-M©) and 
two metallicities (Z = 0.0001 -solid- and Z = 0.0006 - 
dashed-). The boxes corresponding to our selections are 
shown. Both panels show the presence of a sparse se- 
quence above the red HB. While a few of them, which 
are located close to theoretical prediction for the most 
metal-poor stars, fall outside our box, the majority of 
stars, inside the box, look more compatible with slightly 
more metal-rich tracks. Note that this is in agreement 
with the SFH results, which gives for the two galaxies a 
mean metallicity of the order of Z = 0.0004 (Cetus) and 
Z = 0.0006 (Tucana). 

To further investigate the interpretation that the stars 
located inside the box showed in Fig. [3] can be really 
considered the progeny of the stars located along the MS 
blue plume, we compared the star counts ratio between 
the objects located inside the two boxes with theoretical 
predictions obtained ta king into the a ccount the infor- 
mation from the SFH. iMonelh et all (|2010bL their Fig 
18) have already shown a synthetic CMD built from the 
SFH solution of Cetus, where stars with colors and mag- 
nitudes typical of a population of BSSs clearly appear. 
We therefore counted the stars using the same boxes de- 
fined here, and estimated the same ratio using a best-fit 
CMD for both Cetus and Tucana. The derived ratios 
are 0.08±0.01 and 0.08±0.01 for Cetus and Tucana re- 
spectively. Analogous estimates can be made from the 
observed CMD. The number of BSSs in Cetus and Tu- 
cana are 940 and 1214, respectively, while the number of 
objects in the evolved phase are 90 and 123. This gives 
ratios of 0.10±0.01 for Cetus and 0.10±0.01 for Tucana. 
This shows that there is good agreement with the empir- 
ical values. This suggests that the selected stars are the 
evolved counterpart of blue plume objects. Moreover, 
this test strongly suggests that, if we are dealing with 
a pure BSSs population, the analysis performed using 
models calculated for normal MS stars is only marginally 
affecting the results. 

It is also worth noting that the progeny of BSS - in 
particular the more massive ones - can cross the instabil- 
ity strip for radial pulsation during the core He-burning 
stage; in this case they pulsate as anomalous Cepheids 
(jBono et al.lll997t ). When studyi ng the populat i ons o f 
variable stars in the LCID galaxies. iBernard et al.l ()2009l ) 
discovered three anomalous Cepheids in Cetus, and six 
in Tucana. The locations of these variable stars in the 
CMDs of both galaxies are shown in Fig. |3j They ap- 
pear bluer and brighter than the expected bulk of the 
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Fig. 4. — Normalized luminosity function of the candidate BSS population in Cetus (left panel) and Tucana (right panel). The different 
lines refer to the BSS stars within (dotted) and outside (dashed) 1.5r c from the center. 
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Fig. 5. — Ratios of the number of BSSs with respect of HB (asterisks), SGB (full circles), and RGB (open circles) stars, as a function of 
the galactocentric radius for Cetus (left) and Tucana (right). 



BSS progeny. This suggests that anomalous Cepheids 
might have slightly larger masses than the typical BSS. 

3.3. Luminosity function 

The luminosity function of blue plume stars is a pow- 
erful diagnostic to both infer information on their nature 
and, in case of a BSS population, to distinguish between 
the two main formation mechanisms. In fact, if these 
stars belong to an intermediate-age population, we ex- 
pect to find them more centrally concentrated (see next 
section) and also that the brightest (youngest) are prefer- 
entially located to the innermost regions. On the other 
hand, it has been found that in GGCs the BSS lumi- 
nosity function change s as a function of radial distance 
(|Monkman et alj [20061. in the sense that the brightest 
BSS stars are preferentially located in the central regions. 
Since there is no doubt that in GC we only have BSSs, 



with no contamination from intermediate- age stars, this 
occurrence is explained with a combination of two dif- 
ferent formation mechanisms. Collisions, which are effi- 
cient in the inner r egions of GGC, tend to create bluer 
and b righter BSSs iBailvn fc Pinsonn cault 19951 iBailvnl 
1995]). In the outer parts of the GGCs, the BSSs are pre- 
dominantly created from primordial binaries. In dSph 
galaxies the formation of COL-BSS is highly disfavored 
due to the very low stellar density, s o that we do not ex- 
pect to observe a sim ilar behavior (jMapelli et al.l [20071 : 
iMomanv et al.l 120071 ). Therefore, in case a central con- 
centration of the brightest blue plume stars were ob- 
served, this would point to the presence of a genuine 
intermediate-age population. 

Figure [4] shows the BSS candidates luminosity function 
for Cetus (left panel) and Tucana (right panel). The 
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Fig. 6. — Normalized cumulative distributions for the samples of BSSs (black solid line), RGB (red dashed), HB (blue dot-dashed), and 
SGB (green) stars for Cetus (left panel) and Tucana (right). 



plot shows the luminosity function of stars inside (dotted 
line) and outside (dashed line) a distance of 1.5 times 
the core radius, r c . Due to the different number of stars 
and area covered by the two samples, the curves were 
normalized to their area. The shape of the luminosity 
function does not change significantly as a function of 
radius, and similarly we do not find any evidence for an 
obvious shift in the magnitude of the peak. 

We find that the fraction of the brightest BSSs (BS br , 
i.e., within 0.5 mag of the brightest), is also constant 
with radius. For the inner and outer regions we found 
N(BS br )/N(BS) = 0.12 ± 0.02 and 0.12 ± 0.03 in the 
case of Cetus, and N(BS br )/N(BS) = 0.12 ± 0.03 and 
0.13 ± 0.02 for Tucana. Note that the errors are solely 
Poisson statistical errors. 

Thus, we find no evidence for segregation in either 
galaxy, supporting the absence of centrally concentrated 
bright stars, and supporting that the population of BSSs 
in Cetus and Tucana is formed as a result of the evolution 
of binary systems, even in the innermost regions. 

3.4. Radial distribution 

This last result can be further verified studying the ra- 
dial distribution of candidate BSSs compared to that of 
other stellar tracers. It has been found in many GGCs 
that the rad ial profile of BSSs has a central cusp (c.f., 
lBailvnl fl995l). which fall s to a minimum (the 'zone of 
avoidance' IMapelli et al.ll2004f ). and then increases with 
radius. This phenomenon has again been interpreted 
with the different formation mechanisms at work: the 
COL- BSSs are more centrally con centrated than MT- 
BSS (|Bailvn fc Pinsoim cault 1995), which instead dom- 
inate in the outskirts of clusters. If COL-BSSs are not 
present in dSphs, we expect no cen tral concentrat i on, as 
verified bv lMapelli et alJ (|2007t ) and IMapelli et al.l (|2009l ) 
for Draco, Ursa Minor and Sculptor. 

Figure [5] presents the ratios Nbs-HB = Nbs/Nhb, 
NbS-RGB = Nbs/Nrgb, and N bs _sgb — N B s/Nsgb 
for Cetus (left) and Tucana (right), as a function of the 
galactocentric distance expressed in units of core radius. 



We divided the complete catalog into five (Cetus) and six 
(Tucana) radial bins, containing the same total number 
of stars, according to the position of stars in elliptical 
annuli. The points in Figure [5] show the ratio in each 
elliptical radius, plotted as a function of the mean galac- 
tocentric distance of each annulus. 

Figure [5] clearly shows the absence of strong central 
concentrations of BSSs in both galaxies. We tested the 
hypothesis that the ratios of the various samples are con- 
sistent with a flat distribution by comparing the ratios 
to the mean of the observed distribution. In the case of 
Cetus, we estimated non-reduced % 2 of 6.6, 7.3, 9.3 for 
NbS-RGB, Nbssgb, and N bs _hb, respectively. The 
null hypothesis therefore has corresponding probabilities 
of 0.15, 0.12, and 0.05. This suggests that Nbs-RGb and 
Nbssgb are consistent with a flat distribution, while 
NbS-HB is only marginally consistent. In the case of Tu- 
cana, we derive x 2 values of 14.8, 10.6, 12.8 (Nbs-RGB, 
Nbssgb, Nbs-Hb), corresponding to probabilities of 
0.01, 0.06, and 0.03. This means that NbS-RGB and 
NbS-HB are not consistent with a flat distribution, while 
Nbssgb is only marginally consistent. 

Although Nbs_hb has larger uncertainties than the 
other ratios, it is interesting to note that, qualita- 
tively, NbS-HB shows larger fluctuations than the others. 
Moreover, in the range < r c < 2 the trend for the two 
galaxies is opposite. Cetus presents an initial increase 
followed by a drop, while in Tucana an initial decrease 
within the core radius is followed by the increase in the 
ratio of BSSs relative to HB stars. At large distance, 
there is mild evidence of an increase of N B S-HB, espe- 
cially for Tucana where the field covers a large area in 
terms of r c . This is i n agreement with the predictions by 
IMapelli et all (|2006l) . 

Another way to represent the radial distribution of 
different samples is to use the cumulative distributions, 
shown in Figure[6]for the BSS, RGB, SGB, and HB stars. 
A Kolmogorov-Smirnov test gives probabilities of 0.861 
(0.004), 0.48 (0.02), and 0.58 (0.28) for Cetus (Tucana) 
that the BSSs and the RGB/HB/SGB samples are drawn 
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Fig. 7. — The logarithmic frequency of BSSs in Tucana (black 
circle) and Cetus (black diamond) as a function of the absolute 
visual magnitude. We also sh ow the data co r respon ding to other 
LG dwarfs as obtained from Momanv et al. (2007), and Fornax 
(crossed circle, Del Pino et al, 2011, in prep.). 

from the same parent population. 

However, it is interesting to note that, starting from 
3 r c , the BSS cumulative distribution is systematically 
lower, in both galaxies, than the others. This indicates 
higher BSS frequencies at larger distances. 

3.5. Frequency versus My 

iMomanv et al.l (|2007t ) found that a linear correlation 
exists between the logarithmic frequency of BSSs (Fbs, 
normalized to the number of HB stars) and the abso- 
lute visual magnitude of nearby LG dwarfs. Interest- 
ingly enough, they found that this correlation is verified 
over a large range of magnitudes, and that its slope is 
significantly di fferent, i.e., less ste ep, than a similar rela- 
tion found by IPiotto et al] (j2004ft for GGCs. Presently, 
an underlying physical mechanism for this trend has not 
been identified. Based on the selections presented in Fig- 
ure [TJ we estimated that Fbs is equal to 0.05±0.05 and 
0.02±0.04 for Cetus and Tucana, respectively. 

In Figure [3 these values are plotted tog ether with data 
from o ther galaxies given in Figure 2 of Momanv et al. 
(2007) . We adopted the ab solute My magnitude from 
IMcConnachie fc Irwinl (120061) for Cetus (-10.1 mag) and 
from lSaviane et all (|1996l) for Tucana (—9.6 mag). Fig. [7J 
shows that Tucana and Cetus are consistent with the 
relationship derived from oth er LG dwarfs. Using the 
analytic relation for Fbs from IMomanv et al.l (|2007T ) . we 
obtain -0.01±0.08 and 0.03±0.08 for Cetus and Tucana, 
respectively, again, showing very good agreement with 
the observed values. 

If we accept that the Fbs — My relation is signifi- 
cant, then this is an independent indication that the we 
are dealing with a genuine BSS population. If the blue 



plume stars were due to fluctuations in the SFH produc- 
ing younger stars, one would not expect any relationship 
between Fbs and My. On the other hand, if blue plume 
stars are BSSs, they essentialy belong to the same pop- 
ulation generating th e HB stars. As already noted by 
IMomanv et al.1 (|2007l ). the trend is similar in the case of 
globular clusters, but with significantly shallower slope. 
The existence of an anti-correlation with the luminosity 
(the higher the luminosity, the lower the BSS fre- 
quency), and hence with the mass of the cluster, might 
suggest that in more massive globular clusters the mech- 
anisms affecting the binary systems are more efficient 
then in less massive ones, and thi s cause a deficienc y of 
BSS with respect of HB stars (see IPiotto et al.ll2004D . In 
dSph, the shallower slope might indicate that the same 
mechanisms are still at work, but are less efficient due to 
the lower density of dwarf galaxies. 

On the basis of present results, it appears evident that 
if the reliability of Fbs — My relation was confirmed by 
increasing the number of LG dwarfs in the sample, it 
could represent a valuable tool for discriminating stellar 
systems hosting truly young MS stars from those with 
a significant population of old, genuine BSSs. For ex- 
ample, a preliminary analysis of a deep CMD of Fornax 
dSph (Del Pino et al. 2011, in preparation, also shown 
in Figure [7] as a crossed-circle symbol) suggests that the 
Fbs estimate for this galaxy (Fbs > 0.05±0.05) is signif- 
icantly higher than expected from its absolute magnitude 
(Fbs = — 0.23±0.08 ). This is in ag r eemen t with the con- 
clusion presented in iMapelli et a l. (2009), and with the 
fact that Fornax hosts a pop ulation as young as 100 Myr 
(see also lStetson et al.lll998l) . 

4. DYNAMICAL SIMULATIONS 

Further clues about the origin of BSSs can be 
inferred from their dynamical evolution (see, e.g., 
Sigurdsso n et al.l fl994f ) . We thus perform a wide sam- 
ple of dynamical simulations, following the evolution of 
MT-BSS and COL-BSS populations in the potential of 
Cetus and Tucana. 

4.1. Method 

We use the code BEV (Sigurds son fc Phinnevlll995f). in 
the u pgraded version described bv IMapelli et al.T (2004. 
2006). The same code has been used to model BSSs in 
Draco, Ursa Minor (IMapelli et al.|[2007l) . Sculptor and 
Fornax (IMapelli et al.l I2009D . BEV integrates the dy- 
namical evolution of BSSs under the influence of the host- 
galaxy potential, dynamical friction, distant encounters 
with other stars, and three-body encounters (which are 
negligible for dSphs) . The potential of the host galaxy is 
assumed to be time-independent and is represented by a 
multi- mass King model. T he classes of mass are the same 
as m IMapelli et all (120071 ). and the TO mass is assumed 
to be 0.8 Mq (see Figure [2]) for both Cetus and Tucana. 
To calculate the potential, we input the number density 
of stars in the core (n c = 0.03 and = 0.12 stars pc" 3 for 
Cetus and Tucana, respectively) the velocity dispersion 
in the core (we assume a c = 10 km s _1 for both Cetus 
and Tucana^) . and we modify the central adimensional 

15 The observed velocity dispersion, ti f, s , is higher (~ 
17 km s _1 ) for both Cetus IjLewis et al.l 120071 ) and Tucana 
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Fig. 8. — Simulated stellar radial profile (dashed line, red on the web) for Cetus (left-hand panel) and Tucana (right-hand panel), 
compared with observations (filled circles, blue on the web). The number density is g iven in stars per squa re arcmin. Data points for Cetus 
are based on deep Subaru data (Bernard et al. 20ff, in prep.), while are taken from Bernard et al. (2009) for Tucana. 
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Tucana, respectively). Error bars account for Poisson uncertainties. 



potential Wq (denned in iSigurdsson fc Phinnevl Il995f h 
equating the simulated concentration and radial density 
profile with the observed values. Fig. [8] shows the result- 
ing models for Cetus (obtained for Wq = 5 and c = 1.20) 
and for Tucana (obtained for Wq = 1.9 and c = 0.73). 
In the case of Tucana, the simulated profile and the con- 
ce ntration of the model match the observations reported 
in lBernard et a l. (2009). For the concentration of Cetus, 
we adopt c = 1.2, which best matches our data. We note 
that the model for Cetus has a relatively high value of 
Wq, with respect to other dSphs (generally, Wq > 6 for 
globular clusters and Wo < 3 for dSphs). 

In our model, BSSs are generated with a given position, 
velocity, mass and lifetime. The fundamental differences 
between COL-BSSs and MT-BSSs are the following: 

(Fratcrnali et al. 2009), but there is evidence that a significant part 
of it might be due to rotation. 



(i) The distribution of initial positions: the initial po- 
sitions of MT-BSSs are randomly chosen according to a 
probability distribution homogeneous in radial distance 
from the center. This means that MT-BSSs are initially 
distributed according to an isothermal sphere, since MT- 
BSSs are expected to follo w the same distribution a s 
the primordial binaries (see iMapelli et all 12004 [2006h . 
(The spherical nature of the dSphs is not conclusively es- 
tablished, but is chosen as a reasonable starting point.) 
The minimum and the maximum value of the distribu- 
tion of initial radial distances of MT-BSSs, r m j ni mt and 
r msx mt (both defined as three-dimensional quantities), 
have been tuned in order to find the best-fitting simula- 
tion. On the other hand, COL-BSSs are generated only 
inside r c (i.e., r miniC OL = and r maXiC OL = r c in all 
the simulations), the only region were collisions might 
occur, although collisions are highly unlikely in dSphs 
(the probability of a three-body interaction in the core 
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of Cetus and Tucana is > 10 5 times lower than in a mod- 
erately dense globular cluster). The initial distribution 
of COL-BSSs is not necessarily homogeneous in the ra- 
dial distance, as they are not expected to follow the same 
distribution as binaries (in our simulations, we assume a 
constant probability distribution between the center of 
the cluster and r c ). 

(ii) The distribution of initial velocities: the initial ve- 
locities of MT-BSSs are drawn from the equilibrium dis- 
tributions of dSph stars with the same masses as the MT- 
BSSs. C OL-BSSs are assumed to be born with a natal 
kick (see lSigurdsson et al.lll994f ) . We adopt a kick veloc- 
ity, Ukick, = 1 — 2 <r c . For kick velocities higher than 2 o~ c 
most of COL-BSSs are ejected from the dSph, whereas 
for kick velocities lower than 1 <j C) COL-BSSs are almost 
indistinguishable from MT-BSSs from a dynamical point 
of view. 

In most of the runs, the masses of the BSSs are as- 
sumed to be ibbs = 1-3 M©. We also calculated sim- 
ulations with ttibs = 0.9, 1.1 M©, but the results are 
substantially unchanged. Each BSS is evolved for a time 
t, randomly selected from a homogeneous distribution 
between t = and t = tnf e . The param eter t\a e is the 
lifetime of BSSs (see IMapelli et alll2004l) . We calculated 
simulations with t\\f e — 2 and 4 Gyr (4 Gyr is the age 
inferred from the derived SFH). 

4.2. Comparison with observations 

Tables [5] and [3] show the parameters of the simula- 
tions and report also the fraction r\ of COL-BSSs present 
in each run (n = Ncol-bssA^col-bss + ^mt-bss), 
where -/Vcol-bss and -^mt-bss are the number of COL- 
BSSs and of MT-BSSs, respectively). The results of the 
\ 2 analysis are also included. For each run (we made 38 
runs for Cetus and 37 for Tucana) 10000 BSSs have been 
simulated. We show three values of the non- reduced x 2 , 
referring to -/Vbs_rgb (Xrgb' column 8 of Tables [2] and |3]), 
-Nbs_SGB (xIgb> column 9) and iV BS _HB (Xhbj column 
10). Fig. [S] shows the radial distributions of -/Vbs_rgb , 
A'bS-SGB) and -ZVbs_hb obtained for the best-matching 
simulations (run C3 and run T7 for Cetus and Tucana, 
respectively), compared with the observed distributions 
(the same as in Fig. [5]). We stress that only the radial 
distribution of iVess comes from the simulations, and it is 
compared with the observed radial distribution of ^Vrgb , 
TVhb and iVsGB- In fact, the adopted simulation method 
does not allow to evolve these specific stellar types, but 
only stars in a given class of mass (e.g., Mapelli et al. 
2004). We also remind that the underlying potential 
of the dSph is analytic and time-independcntP^I. The 
best-matching simulations reproduce well the observed 
radial distributions of -/V B s_rgb, ^Vbs_sgb, and -/Vbs_hb- 
In fact, the values of the minimum non- reduced % 2 are 

Xrgb ^ 6 > Xsgb ^ 5 and Xsgb ~ 4 for Cetus ( run C3) 
and Xrgb ~ 3 > Xsgb ~ 2 and Xsgb ~ 2 for Tucana (runs 
T5, T7, T14, T15, and T16). 

Fig [TU] shows the behavior of Xrgb (that of Xsgb and 
of Xb_b 1S ver y similar, see Tables [2] and |3]) as a function 

16 A fully N-body simulation, which may account for the dy- 
namical and stellar evolution of all the aforementioned populations, 
would give a much more realistic description of the system, but the 
computational time for running a complete grid of such simulations 
is prohibitive even for star clusters, let alone for dwarf spheroidals. 
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Fig. 10. — Xrgb as a f unc tion of the fraction of COL-BSSs (r/) 
for the simulations of Cetus (squares connected by the dashed line) 
and Tucana (circles connected by the solid line) . For each run, only 
eta was changed, while the other parameters are the same of the 
best models (C3 and T7 for Cetus and Tucana respectively) and 
were held constant. 
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Fig. 11. — The trend with the radial distance of the stellar pop- 
ulation ratios Nggg/Nggg and Nggg /N hhb f° r both Tucana 
and Sculptor (see text for more details). 

of the fraction of COL-BSSs (77). As it was reasonable 
to expect, adding a population of COL-BSSs significantly 
reduces the agreement between data and model: for both 
Cetus and Tucana, only simulations with r\ < 0.2 give an 
acceptable x 2 - 

The dependence of the results on the BSSs mass (in the 
range allowed by the observations) is negligible. This 
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is in agreement wi t h the findings of previous papers 
(|Mapelli et al.ll2007L f2009h for BSSs in Draco, Ursa Mi- 
nor, and Sculptor. The dependence on the lifetime of 
BSSs is negligible for Cetus (where the runs with t\n c — 2 
Gyr are almost equivalent to the runs with fiif e = 4 Gyr) 
and slightly more important for Tucana (where the runs 
with info = 2 Gyr have generally lower values of \ 2 with 
respect to the runs with tuf e = 4 Gyr, but the difference 
is not statistically significant). Given the range of masses 
expected for the BSS, suggested by the comparison with 
stellar tracks, we assume 4 Gyr as fiducial value. 

In the case of Cetus, for the best-matching runs r m - m ~ 
0.0 and r max ~ 5r c , indicating that primordial binaries 
were initially distributed ranging from the center of the 
dSph out to (at least) 5 r c . The tidal radius of Cetus (de- 
rived from our data) is 24', i.e., ~ 15 r c , much larger than 
r max . However, this difference likely does not mean that 
BSSs do not form beyond 5r c , but is a consequence of 
the fact that we do not have data beyond » 4 r c (further- 
more, the coverage of Cetus is not complete even inside 
4r c ). In the case of Tucana, for which data are available 
almost up to the tidal radius, the best-matching value of 
r mm is equal to the tidal radius. Interestingly, the best- 
matching values of r m ; n for Tucana are all larger than 
one core radius. 

In summary, the results of dynamical simulations show 
that it is possible to explain the observed distribution of 
BSS candidates, in both Cetus and Tucana, with a MT- 
BSS population. Any population of COL-BSSs more nu- 
merous than the 10% of the total number of BSSs reduces 
the agreement between data and simulations. This result 
is not surprising, as basic dynamical calculations tend to 
exclude the existence of COL-BSSs in dSphs, but it is 
important to point out that the simulations can almost 
rule out the unphysical scenario where > 20% of BSSs 
are COL-BSSs in both Cetus and Tucana. 

5. COMPARISON WITH THE SCULPTOR DSPH 

The results discussed in the previous sections provide 
evidence that the two isolated dwarf galaxies Cetus and 
Tucana host a significant population of BSSs, whose 
properties strongly suggest that they are the product of 
mass exchange in binary systems. We did not find any 
evidence supporting the hypothesis that part of the ob- 
jects populating the blue plume are intermediate-age or 
young MS stars. On the basis of these findings, it appears 
that b oth Cetus and Tuca na are similar t o Ursa Minor, 
Draco (|Mapelli et al.ll2007f ). and Sculptor (jMapelli et al.l 
2009). In the following, we first focus on the compari- 
son of Tucana and Sculptor, two dwarfs which show close 
similarities in terms of their SFH and present-day stellar 
content. 

Sculptor is a well established case of a dSph galaxy 
hostin g two different old populations. Maicwsk i et al.1 
(1999) fir st noted evidence of two RGB bumps in 
Sculptor. iTolstov et al.l (|2004l) presented a photometric 
and spectroscopic (235 RGB stars) investigation of this 
galaxy, showing the presence of strong radial gradients 
and two stellar populations with different kinematic and 
chemical properties. Moreover, the radial distributions of 
the red and blue HB components are correlated with the 
properties of the RGB stars, suggesting that red stars are 
more centrally concentrated, more metal-rich, and have a 
higher velocity dispersion, while the blue component has 
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Fig. 12. — Ratios as a function of galactocentric distance for 
different samples of stars. BSSs (left) and HB stars (right) are 
compared to red (assterisks, dashed line) and blue (open symbols 
and solid line) RGB stars for Tucana (top) and Cetus (bottom). 

a broader distribution, lower metallicity, and a smaller 
v elocity dispersio n. 

iBernard et al.l ([2008D showed the first evidence of stel- 
lar radial gradients in Tucana; after dividing both the HB 
and the RGB into red and blue components, the num- 
ber of stars belonging to the red components decreases 
faster than the corresponding blue sequence with increas- 
ing galaxy radius. Moreover, t he analysis of the pr oper- 
ties o f RR Lyrae variables s tars ([Bernard et al.ll2~009D , the 
SFH (iMonelli et al.ll2010cD . and the double RGB bump 
(jMonelli et al.ll2010al) strongly support the coexistence of 
two old (> 10 Gyr) populations in Tucana, with slightly 
different ages, metallicities, and spatial distributions. 

Figure [TT] compares the ratio of the number of BSSs 
versus red (asterisks) and blue (squares) HB stars as 
a function of galactocentric radius for the two galax- 
ie s. The values for S culptor are taken from Table 2 
of IMapelli et al.l ((2009). Interestingly, in both cases the 
BSSs are more centrally concentrated than the blue com- 
ponent of the HB, and less concentrated than the red 
component. We note also that, for each galaxy, the two 
curves cross at a similar galactocentric distance, « 1 — 1.6 
r c . In the case of Tucana, this trend is confirmed also for 
the BSSs and the HB stars when compared to red and 
blue RGB stars (Figure[T2J top left and top right, respec- 
tively). Therefore, if we compare the global Nbss/Nrb 
ratio in Figure [5] with Figure ITT1 and Figure [TJ] (top pan- 
els), it can be inferred that the radial trends of the BSSs 
are indeed driven by the fast decrease of the red HB 
stars in the central regions, associated with the younger, 
more metal- rich component, while the increase at larger 
radii is driven by the more homogeneously distributed 
bluest stars. This means that the peculiar radial trend 
of the candidate BSSs when they are normalized to the 
HB stars is not due to a peculiar distribution of these 
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TABLE 2 

Parameters and values of x 2 for the simulations of Cetus. 



CelSG a 


/ b 
r min, MT/ r c 


/ b 
r max, MT/ r c 


m BSS/ 1V1 Q 


life/ *-*y r 


1 'kick/ a 'c 


'/ 


v 2 s 
argb 


v 2 S 
ASGB 


v 2 
AHB 


Cl 


0.0 


5.0 


0.9 


4 


_ 





7.7 


6.5 


4.7 


02 


0.0 


5.0 


1.1 


1 


_ 





9.5 


8.1 


5.8 


C3 


0.0 


5.0 


1.3 


4 


_ 





5.8 


4.9 


3.5 


C4 


0.1 


5.0 


1.3 


1 


_ 





8.3 


7.0 


5.0 


C5 


0.2 


5.0 


1.3 


4 


_ 





7.3 


6.2 


4.5 


C6 


0.3 


5.0 


1.3 


1 







10 


8.7 


6.2 


C7 


0.1 


5.0 


1.3 


1 


_ 





11 


12 


8.5 


08 


0.5 


5.0 


1.3 


1 


_ 





19 


16 


12 


C9 


1.0 


5.0 


1.3 


1 


_ 





51 


12 


31 


CIO 


0.0 


5.0 


0.9 


2 


_ 





6.7 


5.6 


4.1 


Cll 


0.0 


5.0 


1.1 


2 







7.2 


6.1 


1.4 


C12 


0.0 


5.0 


1.3 


2 


_ 





9.8 


8.3 


5.9 


C13 


0.1 


5.0 


1.3 


2 







9.0 


7.6 


5.4 


C14 


0.2 


5.0 


1.3 


2 







7.0 


5.9 


4.2 


C15 


0.3 


5.0 


1.3 


2 


_ 





7.5 


6.3 


4.6 


C16 


0.1 


5.0 


1.3 


2 


- 





12 


10 


7.5 


C17 


0.5 


5.0 


1.3 


2 


_ 





20 


17 


12 


C18 


1.0 


5.0 


1.3 


2 


_ 





50 


11 


30 


C19 


0.0 


4.0 


1.3 


1 


- 





10 


8.7 


6.2 


C20 


0.2 


4.0 


1.3 


1 


_ 





9.4 


7.8 


5.5 


C21 


0.5 


4.0 


1.3 


4 


_ 





7.9 


6.7 


4.8 


C22 


1.0 


4.0 


1.3 


4 







11 


34 


25 


C23 


0.0 


6.0 


1.3 


1 


_ 





13 


11 


7.7 


C24 


0.2 


6.0 


1.3 


1 


_ 





16 


13 


9.1 


C25 


0.5 


6.0 


1.3 


4 


_ 





26 


22 


16 


C26 


1.0 


6.0 


1.3 


1 


- 





77 


64 


16 


C27 


0.0 


5.0 


1.3 


1 


1 


0.1 


7.2 


6.0 


4.3 


C28 


0.0 


5.0 


1.3 


1 


1 


0.2 


13 


11 


7.8 


C29 


0.0 


5.0 


1.3 


1 


1 


0.3 


20 


17 


12 


C30 


0.0 


5.0 


1.3 


1 


1 


0.4 


28 


23 


16 


C31 


0.0 


5.0 


1.3 


1 


1 


0.5 


43 


36 


26 


C32 


0.0 


5.0 


1.3 


1 


1 


0.6 


56 


47 


33 


C33 


0.0 


5.0 


1.3 


1 


1 


0.7 


76 


63 


15 


C34 


0.0 


5.0 


1.3 


1 


1 


0.8 


94 


78 


56 


C35 


0.0 


5.0 


1.3 


1 


1 


0.9 


110 


91 


65 


C36 






1.3 


1 


1 


1.0 


122 


101 


72 


C37 


0.0 


5.0 


1.3 


1 


2 


0.1 


9.4 


7.9 


5.7 


C38 


0.0 


5.0 


1.3 


4 


2 


0.2 


11 


12 


8.3 



a Identifier of the run (for each run, 10000 BSSs have been simulated). 

fmiii, mt and r maXi mt are the minimum and the maximum three-dimensional radius within which MT-BSSs are generated in the 
simulations. The minimum and the maximum three-dimensional radius within which COL-BSSs arc generated in the simulations are 
not listed in the Table, because they are the same in all the simulations (i.e., r m in, col — and r maXi col — r c , respectively). 
c Mass of a simulated BSS. In the simulations, mess is the same for MT-BSSs and COL-BSSs. 

Lifetime of a simulated BSS. In the simulations, tn{ e is the same for MT-BSSs and COL-BSSs. 
c Kick velocity for COL-BSSs. 

^ Fraction of simulated COL-BSSs with respect to the total population of BSSs. 

S Xrgbj Xsgb and Xhb aro tnc non-reduced x 2 of the simulated JVbs_RGB, -WbS-SGB and w bs_hb- 



objects inside the galaxy but only to the existence of 
a significant metallicity gradient in Tucana, a gradient 
that is well traced by both the red and blue HB stellar 
populations. 

The similarities between these two galaxies are par- 
ticularly intriguing because of the remarkably different 
environmental conditions: Tucana spent most of its life 
in isolation, while Sculptor is a close satellite of the Milky 
Way. If all the BSSs are descendants of primordial bina- 
ries, this suggests that the different environments of the 
two galaxies did not alter their binary distributions. 

In contrast, Figure [T2l (bottom panels) shows that the 
properties of different stellar samples in Cetus appear 
more homogeneous. This would be in agreement with 
other findings, as the properties of the RR Lyrae vari- 
able stars and of the RGB bump, which indicate a single 
dominant population characterized by a more uniform 
distribution of ages and metallicities. 

Interesting features appear when comparing the distri- 



bution of the BSSs with those of the red and blue stars 
of both the RGB and HB. We defined these samples as 
follows. The region selecting the RGB stars was sim- 
ply split in two, similarly to what has been done f or the 
analysis of the RGB bump in lMonelli et all (|2010aD . The 
red and blue HB are defined considering only the stars 
redder than m,FA75W — mFsuw > 0.99 and bluer than 
mF475W — mF8i4W < 0.46 in the HB region, therefore 
mostly avoiding the RR Lyrae instability strip. Figure 
1131 compares the distributions of these four samples with 
the BSSs for both Tucana and Cetus. In the case of Ce- 
tus, we find that the BSSs are significantly less concen- 
trated than the red HB and more concentrated than the 
blue HB, while they more closely follow the distributions 
of both of the RGB components. In the case of Tucana, 
the BSSs closely follow the distributions of both blue 
components (possibly more concentrated in the central 
regions), while they are significantly less concentrated 
than the red sequences. 



The BSS population in Cetus and Tucana 
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TABLE 3 

Parameters and values of \ 2 for the simulations of Tucana. 



Case a 


' mm, MT / ' c 


T' -\ Tm It* ^ 

' max, MT/ ' c 


YY1 n n r-i 1 A// ' r-* 

'''■BSS/ ivl Q 


c hfe/ 


/ e 
^kick/ °"c 


n f 
'/ 


v 2 s 

ARGB 


v 2 s 
ASGB 


v 2 > 
AHB 


Tl 


0.0 


5.4 


1.3 


4 





_ 


101 


78 


60 


T2 


0.5 


5.4 


1.3 


4 





_ 


51 


40 


30 


T3 


1.0 


5.4 


1.3 


1 





_ 


21 


17 


12 


T4 


1.5 


5.4 


1.3 


1 





_ 


5.7 


4.6 


3.4 


T5 


1.7 


5.4 


0.9 


1 





_ 


2.6 


2.1 


1.6 


T6 


1.7 


5.4 


1.1 


4 





_ 


3.2 


2.6 


1.9 


T7 


1.7 


5.4 


1.3 


4 








2.4 


1.9 


1.5 


T8 


2.0 


5.4 


1.3 


4 





_ 


7.7 


6.1 


4.8 


T9 


2.5 


5.4 


1.3 


1 





— 


18 


14 


11 


T10 


0.0 


5.4 


1.3 


2 





_ 


106 


83 


63 


Til 


0.5 


5.4 


1.3 


2 





_ 


55 


43 


33 


T12 


1.0 


5.4 


1.3 


2 





_ 


23 


18 


11 


T13 


1.5 


5.4 


1.3 


2 





_ 


2.5 


2.1 


1.5 


T14 


1.7 


5.4 


0.9 


2 





_ 


2.3 


1.8 


1.4 


T15 


1.7 


5.4 


1.1 


2 







1.9 


1.5 


1.1 


T16 


1.7 


5.4 


1.3 


2 





_ 


1.6 


1.3 


0.95 


T17 


2.0 


5.4 


1.3 


2 





_ 


4.8 


3.7 


2.9 


T18 


2.5 


5.4 


1.3 


2 





_ 


16 


13 


9.8 


T19 


0.0 


4.0 


1.3 


1 





_ 


139 


109 


83 


T20 


0.5 


4.0 


1.3 


1 





_ 


83 


65 


49 


T21 


1.0 


4.0 


1.3 


4 







11 


32 


21 


T22 


1.5 


4.0 


1.3 


1 





_ 


15 


12 


8.6 


T23 


1.7 


4.0 


1.3 


1 





_ 


5.5 


4.5 


3.3 


T24 


2.0 


4.0 


1.3 


1 





_ 


6.3 


1.9 


3.8 


T25 


2.5 


4.0 


1.3 


4 





_ 


11 


8.7 


6.7 


T26 


1.7 


5.4 


1.3 


1 


1 


0.1 


7.8 


6.0 


4.5 


T27 


1.7 


5.4 


1.3 


1 


1 


0.2 


21 


18 


14 


T28 


1.7 


5.4 


1.3 


1 


1 


0.3 


45 


34 


26 


T29 


1.7 


5.4 


1.3 


1 


1 


0.4 


98 


76 


58 


T30 


1.7 


5.4 


1.3 


1 


1 


0.5 


140 


109 


84 


T31 


1.7 


5.4 


1.3 


1 


1 


0.6 


199 


156 


120 


T32 


1.7 


5.4 


1.3 


1 


1 


0.7 


247 


195 


149 


T33 


1.7 


5.4 


1.3 


1 


1 


0.8 


297 


236 


180 


T34 


1.7 


5.4 


1.3 


1 


1 


0.9 


330 


262 


200 


T35 






1.3 


1 


1 


1.0 


377 


301 


229 


T36 


1.7 


5.4 


1.3 


1 


2 


0.1 


5.8 


1.5 


3.3 


T37 


1.7 


5.4 


1.3 


1 


2 


0.2 


24 


18 


15 



a Identifier of the run (for each run, 10000 BSSs have been simulated). 

^min, MT and r maXi mt are the minimum and the maximum three-dimensional radius within which MT-BSSs are generated in the 
simulations. The minimum and the maximum three-dimensional radius within which COL-BSSs arc generated in the simulations are 
not listed in the Table, because they arc the same in all the simulations (i.e., r m i n ,cOL — and r maXl COL — r c , respectively). 
c Mass of a simulated BSS. In the simulations, mess is the same for MT-BSSs and COL-BSSs. 

Lifetime of a simulated BSS. In the simulations, infe is the same for MT-BSSs and COL-BSSs. 
c Kick velocity for COL-BSSs. 

Fraction of simulated COL-BSSs with respect to the total population of BSSs. 
S Xrgb> Xsgb and Xhb aro tnc non-reduced \ 2 of the simulated Af B s_RGB, -Wbs.SGB and AT B S_HB. 



We have demonstrated that the populations of BSS 
candidates in Tucana and Cetus are consistent with all 
predictions of creation from MT-BSSs. In fact, it would 
be unprecedented to have a truncation in star formation 
without an associated population of MT-BSSs. Unfortu- 
nately, our limited knowledge of the evolution of binary 
stellar systems prevents us from making accurate predic- 
tions of the expected numbers of MT-BSSs. In sum, we 
expect to sec MT-BSSs in dSph galaxies, and we detect 
them with consistent properties. The remaining question 
is whether some fraction of the BSS candidates could still 
be due to newly formed MS stars. The inconsistency with 
the age-metallicity relationship could be due to infall of 
a gas cloud with a low metallicity coincident with that 
of the main population which then produces subsequent 
star formation. However, the truncation at luminosities 
corresponding to twice the TO mass implies a second 
unlikely coincidence if some of the BSS candidates were 
true MS stars. While very difficult to rule out such a sce- 



nario completely, it is clear that newly formed MS stars 
are not necessary to explain all of the observations. 

6. SUMMARY AND CONCLUSION 

In this paper, we have presented the analysis of can- 
didate BSSs in two isolated dSphs in the periphery of 
the Local Group, Cetus and Tucana. Deep HST/ACS 
data allowed us to identify 940 and 1214 BSS candidates 
in Cetus and Tucana, respectively. Different indicators 
suggest that these are true BSSs; and, thus, there are no 
young main sequence stars in either galaxy. The first ev- 
idence comes from th e SFH we derived for both objects 
(|Monelli et al.l 2010b. c). Not allowing for the presence 
of BSSs, then a small population of relatively young (3 
to 5 Gyr old) stars is detected, with metallicities simi- 
lar to the oldest most metal-poor stars, and thus, much 
lower than expected from the age-metallicity relation of 
the dominant population. The comparison with stellar 
tracks presented in this paper suggests that the BSSs 
have typical masses > IMq, and smaller than ~ 1.5M©, 



14 



Monelli et al. 




Fig. 13. — Normalized cumulative distributions for the samples of BSSs (black solid line) compared to red and blue RGB stars (red and 
blue solid lines), and red and blue HB stars (red and blue dashed). In the Case of Cetus (left panel), the BSSs are more concentrated than 
the blue HB stars, and less concentrated than the red HB stars. In the case of Tucana (right panel), the BSSs are less concentrated than 
both the red RGB and HB samples, while they are marginally more concentrated than the blue components. 



in agreement with the fact that the maximum expected 
mass of a BSS is twice the mass presently evolving at 
the TO (~ 0.8 M Q ). Moreover, the limited spread in 
metallicity found for the BSS populations, perfectly in 
agreement with that of turn-off and subgiant stars, sug- 
gests that the use of stellar evolution models for single 
stars does not introduce strong biases in the analysis. 

The analysis of the luminosity functions and radial dis- 
tributions support the BSS interpretation. On one hand, 
we showed that the luminosity function does not change 
as a function of radius, and in particular the bright- 
est (youngest) stars do not appear more centrally con- 
centrated, as one would expect from a young main se- 
quence. On the other hand, the radial profiles do not 
show any central peak, as observed in globular clusters 
due to COL-BSS stars or as it would be in the case of a 
young component. 

Thus, the present analysis also suggests that the BSSs 
in Cetus and Tucana formed from primordial binaries 
and, as in the case of other dSphs, COL-BSSs are un- 
likely to form. This conclusion is strongly supported by 
the dynamical simulations we performed, taking into ac- 
count the properties of the dwarfs under scrutiny. These 
simulations indicate that the observed BSS distributions 
are well reproduced when only MT-BSSs are taken into 
account, while the fit worsens when even a small num- 
ber of COL-BSSs (~ 10%) is included in the dynamical 
simulation. 

We can therefore safely conclude that it is highly likely 
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